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ABSTRACT 

Approximate equations which describe the behavior of cosmic rays 
in the interplanetary medium under suitable conditions are used to make 
comparisons between observations and theoretical predictions of radial 
gradients and radial anisotropies. In the high energy region there appear 
to be no inconsistencies between theory and observations. In the low 
energy region it is shown that theoretical predictions of the radial 
anisotropy expected from large radial gradients of the intensity are not 

inconsistent with observed radial anisotropies. However, in the latter 

i ' 

case there are other inconsistencies , which suggest that some aspects 

i 

of the observations or of the theory (or both) are unsatisfactory. 

1 . Introduction 

The modulation of cosmic ray particles in the interplanetary 
medium can be discussed in terms of the quasi-steady, spherically- 
symmetric model developed by Parker (1965, 1966), Gleeson and Axford 
(1967, 1968a), and Jokioii and Parker (1967). In this paper observa- 
tions of radial gradients and radial anisotropies of both high energy 
and low energy cosmic rays ate compared with the predictions of the 
approximate equations derived for this model by Gleeson and Axford 
(1968b, c) and by Fisk and Axford (1969). 

It has been shown by Gleeson and Axford (1967) that the cosmic ray 
number density U(r,T) and streaming (or radial current density) S(r,T), 

A* ' 

per unit interval of kinetic energy T, satisfy the equations: 
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is the Comp ton-Get ting factor (Gleexon and Axford., 1968a). Here r is the 
heliocentric distance, K (r,T) is the diffusion coefficient, V(r) is the 
solar wind speed, and a(T) ■ (T + 2T 0 )/(T + T 0 ) where T 0 is the rest 
energy of a particle. When V is a constant, S can also be expressed as: 

r a au 

s - *2 ta (vu - K a? > w 


Analytic solutions of these equations have been found for certain 
simple forms of K (r,T), with a assumed constant, and (usually) with the 
unmodulated cosmic ray spectrum a power law in kinetic energy (see Fisk 
and Axford, 1969). Numerical solutions are available which do not suffer 
from these restrictions, although they must inevitably involve arbitrary 
assumptions about the variation of K with r and the form of the unmodulated 
spectrum (Fisk, 1969). Unfortunately, while these solutions are valuable 
in a qualitative sense, they are not directly useful in interpreting obser- 
vations of cosmic rays in the vicinity of the earth. It is more convenient 
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for this purpose to use simpler, more manageable approximate equations 
which are valid in some limited energy range and over some part, if 
not all, of the modulating region. Useful approximations are possible 

f* 

in cases where the dimens ionaless number ft 2 V r/K (with the tilde 
denoting 'characteristic value') is either sufficiently small or 
sufficiently large (Gleeson and Axford, 1968b, c; Fisk and Axford, 1969). 
The analytic and numerical solutions described above are useful in that 
they permit the range of validity of these approximate equations to be 
estimated. 

In this paper we discuss these approximate equations further, 
emphasising in particular their usefulness for determining the behav- 
ior of the radial anisotropy. In Section 2 we discuss the approximate 
equation for the number density which is valid when ft 1, a condition 
that should be satisfied in the vicinity of the earth at energies above 
a few hundred MeV/nucleon. We find that solutions to this equation can 
be used to determine corresponding approximations to the streaming and 
hence the radial anisotropy only to limited accuracy. In Section 3 we 
discuss the approximate equations relating the anisotropy to the spectrum 
when ft is large, which is likely to be the case at energies below 50-75 
MeV/nucleon. Assuming several possible forms for the low energy proton 
spectrum, we use these equations to determine the corresponding aniso- 
tropies, which are then compared with the anisotropies observed at low 
energies by Rao, et al. (1967). We find that there is an acceptable 
spectrum for which the observed and predicted anisotropies agree. 
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Finally, in Section 4 we review earlier studies of the behavior of the 
anisotropy . 

2. Approximations for Intermediate and High Energies 

It has been argued by Gleeson and Axford (1967, 1968c) and by Fisk 
and Axford (1969) that if only galactic cosmic rays are considered, 
the streaming S can be neglected on the left side of Equation (2) when 
|(l/2) C (C - l)ft|« 1. In this case the number density satisfies the 
* force- field 3 equation: 

CVU =“ K ~ (5) 

or 

It is believed that this approximation j.s valid locally (i.e. in the 

H 

vicinity of the orbit of the earth) for particles with energies exceed- 
ing a few hundred MeV/nucleon, and beyond the orbit of the earth for 
particles with even lower energies. In the energy range for which this 
approximation is valid, Equation (5) relates the radial gradient of the 
number density directly to the spectrum of the cosmic ray particles and 
to the diffusion coefficient observed locally. It should be noted that 
the factor C in the expression for the gradient has an Important effect 
at energies below about 1 GeV/nucleon, and makes the gradient in this 
energy range much smaller than would be expected on the basis of the 
simple convection-diffusion model (which corresponds formally to C * 1). 

An approximate expression for the radial anisotropy, § * 3S/vU (where 
v is particle speed), can be obtained by substituting Equation (5) into 


Equation (A); thus 
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and hence 
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In fact, this is only an order- of -magnitude result and should generally 
be expected to hold to within a factor ~ 2. The reason for this numerical 
uncertainty can be seen by comparing Equation (6) with Equation (1). It 
is evident that Equations (6) and (7) can be expected to be good approxi* 
mations only if dS/dr « S/r; however, it can be shown from analytic 
solutions that ds/dr is in general of the same order as S/r when ft < 1, 
and hence an error of a factor ^2 can be expected. For the very simple 
solution in which k * K 0 r, Equation (6) yields the correct asymptotic 
result in the limit 3 « 1, but when Kwr 1 *, the error appears to be a 
factor 2/(3 - b) . 

To determine § to lowest order in & it is clear from Equation (7) 
that it is sufficient to use the unmodulated spectrum for C and U. Thus 
if we take U 06 T*^ at high energies, then C » [l + Of(p. - 1) /3 ]. Further- 
more, there is some observational evidence which suggests that for high 
energies the diffusion coefficient has the form K » 0PfC.i(r), where 0 « v/c, 
c is the soeed of light, and P is the particle rigidity (Gloeckler and 
Jokipii, 1966). Noting that 0P 06 orT we find that with these forms for U, 

G, and k. 


* 
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According to O'Gallagher (1967) the radial gradient of protons with 
T ~ 5 GeV in the heliocentric distance range 1 • 1.57 A. U. wee 
(r/U) (dU/dr) ** +97. daring 1965. On taking V ■ 400 km. sec."* - 0.0014v, 

0t - 1.17, and - 2.65, we find from Equation (8) that § *■ +0.027.. Thia 
small positive anisotropy at high energies is in agreement with earlier 
predictions by Gleeson and Axford (1968b, c). It is sufficient to chift 
the direction of maximum intensity associated with the normal diurnal 
variation from 90° east to about 88° east of the sun, in accordance with 
observations described by McCracken and Rao (1966). 

In Figure 1 we have plotted the radial gradient of the Intensity 
determined by Equation (5), and the radial anistropy determined by 
Equation (7), using the spectrum for protons given Gloeckler and Jokipii 

* 

(1967) and assuming that K « |3p. ! 

i 

i- 

| 

3 . Approximations for Low Energies j 

i. 

It seems well-established that in the vicinity of the earth the 
diffusion coefficient decreases with decreasing energy at least down i 

i 

to energies of the order of a few hundred MeV/nucleon (e.g. Ormes and [ 

Webber, 1968; Sari and Ness, 1969). If this behavior continues we [ 

anticipate that locally ft attains quite large values (i.e. ® £ 3) at 1 

energies below about 50 - 75 MeV/nucleon. * 

Fisk and Axford (1969) have shown that two distinct sets of approxi- j 

't 

mate equations are possible in the limit of large ft: one is applicable if 
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the magnitudes of the radial gradients of the number density and stream- 
ing are of order ft , and the other, for substantially smaller gradients, 
Unfortunately, it is not possible to decide at present which of the two 
sets should be used since radial gradients of the intensity with both 
large and small magnitudes have been observed at low energies (e.g. 
O'Gallagher, 1967; Anderson, 1968; Krimigis and Venkatesan, 1968), We 
will discuss the local particle behavior in termp of both sets of pos- 
sible approximate equations valid for large ft . 

Fisk and Axford (1969) showed that if the gradients at low energies 
are as large as observations by O'Gallagher (1967) indicate, then the 
local behavior of galactic cosmic rays at energies below 50 - 75 HeV/ 
nucleon can be described by an approximate equation which relates the 
streaming directly to the local spectrum: 

S ~ “ 1 JL (OTU) (9) 

3 dT 


This equation can be obtained directly from Equation (1) by dropping 
the term 2s/r and integrating, assuming V to be constant. The inte- 
gration produces an additional function of T, which we interpret as 
representing a source term (i.e. solar particles), and accordingly 
neglect. (The term is indeed absent in our analytic solutions (Fisk 
and Axford, 1969).) In terms? of the radial anisotropy § and the mean 
differential intensity jo ■ vU/4rr, Equation (9) can be written 



In jo 
In T 
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S 


( 10 ) 
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In this case the number density of galactic cosmic rays satisfies the 
simple convection-diffusion equation: 

VU (11) 

bt 


provided the spectrum does not have an extremely large slope, that 
is, provided 

I" ’ i I “ttH* + l l ~ 1 <12) 

If the radial gradients at low energies are small (as suggested by 
Anderson (1968) and by Krimlgls and Vehkatesan (1968), for example), 
then if ft is large, it is permissible to neglect the term K dU/dr on 
the right side of Equations (2) and (4), thus yielding: 


and 


S CVU 
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In terms of the radial anisotropy and the mean differential intensity, 
Equation (13) can be written: 


« - - 1 < • fiH 2 • 2) (15) 

On eliminating S between Equations (13) and (14) we obtain a simple 
expression for the radial gradient: 


JL £io =- 2 (<* 8 In jo 
jo or 3 a in T 
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This approximation corresponds to the cese in which the scattering la so 
effective that it keeps the cosmic ray particles "frozen" in the expand- 
ing solar wind, and hence they undergo severe deceleration. Equation (15) 
accordingly corresponds to the anisotropy associated with a bulk speed V, 
allowing for the Compton-Getting effect (Gleeson and Axford, 1968a, 

Forman, 1969), while Equation (16) corresponds to pure adiabatic expan- 
sion. These approximate equations can be used to treat galactic cosmic 
ra>s (Fisk and Axford, 1969), and also presumably solar cosmic rays (e.g. 
Forman, 1968a) provided | (r/j 0 ) dj 0 /dr| « ft. 

Having derived these two sets of approximate equations for the case 
of large ft, we are in a position to discuss the implications of obser- 
vations of radial gradients and anisotropies at low energies. In 
Figure 2 we have plotted a spectrum which provides a reasonable fit 
to the available observations of low energy proton intensities in 1966. 
This spectrum could presumably result from (a) galactic protons for 
which the gradients are large, (b) a combination of galactic protons, 
for which the gradients are small, together with solar protons, or (c) 
a combination of galactic protons, for which the gradients are large, and 
solar protons. To determine which, if any, of these three possible cases 
is acceptable (assuming that locally ft is large at low energies), we can 
compare the anisotropy determined for each of the cases by Equations (10) 
and (15) with the low energy anisotropy observed by Rao, et al. (1967). 
The particles observed by Rao, et al. in the vicinity of the earth 

(r * 0.8 - 1.1 A. U.) in 1966 are presumably mainly protons, and have 
essentially radial mean anisotropies of § « 0.19 + 0.067# in the energy 
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range (7.5 • 45 MeV), and of § * 0.2 + 0.2% in the range (45 - 90 MeV). 

The approximate equations are unlikely to be valid at energies above 
50 • 75 MeV, and hence a comparison between predicted and observed 
anisotropies in the (45 - 90 MeV) energy range is probably not mean- 
ingful. We shall, however, give the predicted anisotropies in this 
upper energy range for completeness. 

Case (a): Galactic protons for which the gradients are large . 

The radial anisotropy determined by Equation (10) for the form of 
the spectrum shown in Figure 2 is plotted in Figure 3. To determine 
the predicted mean anisotropies, we average the anisotropy shown in 
Figure 3 over each of the two observed energy ranges using the inten- 
sity as a weighting function: 

T3 Tg 

5 - j § Jo «/ f Jo dT (17) 

r Ti 

where Tg and Ti are the upper and lower limits, respectively, of a given 
energy range. For this case the predicted and observed anisotropies agree 
in the (7.5 - 45 MeV) energy range where the predicted mean anisotropy is 
found to be f ■ 0.25%. In the (45 - 90 MeV) energy range the predicted 
mean anisotropy of ^ ■ -1.14% is much larger in magnitude than the observed 
anisotropy, but, as we indicated above, close agreement should not be 
expected in this energy range. In Figure 4 we have plotted the quantity T| 
(Equation (12)) for this case. As is evident in this figure, T) is gen- 
erally of the order unity or less; hence the number density corresponding to 
this form of the spectrum should satisfy a simple convection- diffusion 
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quation at low energies, and V/K can be estimated from the radial 
gradient . 

Case (b) ; Galactic protons for which the gradients are small, together 
with solar protons. 

The radial anisotropy determined by Equation (15) for this case is 
also plotted in Figure 3. Here the predicted mean anisotropies are 
found to be f ■ 2.30% in the (7.5 - 45 MeV) energy range, and f « 0.0001% 
in the (45 * 90 MeV) range. The predicted and observed anisotropies do 
not agree in the lower energy range, and such agreement as exists in 
the uoper energy range is probably fortuitous. In Figure 5 we have 
plotted the radial gradient of the intensity determined by Equation (16) 
for this form of the spectrum. The pronounced negative gradient at energies 
below 65 MeV shown in Figure 5, and the large positive anisotropy shown in 
Figure 3 at these energies presumably indicate that this form of the 
proton spectrum could only be realistic if the protons at low en:: 0 ies 
were predominantly of solar origin. Since the observed and predicted 
anisotropies do not agree in the lower energy range, this is apparently not 
the case. It should be noted that at very low energies (< 12 MeV) the 
magnitude of the radial gradient shown in Figure 5 is large, and hence it 
is not certain whether the condition required for Equations (15) ord (16) 
to be valid, viz. (1/ft) | (r/ j 0 ) (d j G /dr) J 1, is satisfied. We assume, 
however, that these equations indicate at least the general features of 
the particle behavior predicted by this form of the spectrum. 
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Case (c): Galactic protons for which the gradients are large, together 

with solar protons. 

For this case we assume that the spectrum shown in Figure 2 can be 
decomposed into solar and galactic proton spectra, each of which is 
also shown in Figure 2. The solar proton anisotropy is determined by 
Equation (15), and the galactic proton anisotropy by Equation (10). 

The predicted anisotropy shown in Figure 3 is then obtained by adding 
each of the anisotropies, weighting them according to the relative 
number of particles of each species present at a given energy. The 
predicted mean anisotropies of §" « 1.277. in the (7.5 - 45 MeV) energy 
range, and of X m -1.12% in the (45 - 90 MeV) range are much larger 
in magnitude than the corresponding observed anisotropies. 

According to the above discussion, of the three cases considered, 
only the anisotropy corresponding to the form of the low energy proton 
spectrum given in case (a) (galactic protons for which the gradients 
are large) fits the observed anisotropy. Of course, the predicted 
anisotropies are sensitive to the shape of the curve used for the 
spectrum, which, since there are only a few data points available in 
the energy range considered (7.5 - 90 MeV), is not well established. 
However, with reasonable variations of the spectrum it is unlikely 
that the anisotropies in cases (b) and (c) would agree with the ob- 
served anisotropy. The mean anisotropy in case (a), even with some 
variation of the spectrum, should be small and positive in the (7.5 - 
45 MeV) energy range, in agreement with the observations. 
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We conclude, therefore, that (1) a low energy proton spectrum 
which is the result primarily of galactic protons and which resembles 
the spectrum shown in Figure 2, (li) a radial gradient of the inten- 
sity of galactic orotons which locally has a large magnitude at low 
energies (i.e. a gradient as large as that observed by O'Gallagher 
(1967)), and (iii) the small low energy anisotropy observed by Rao, 
et al. (1967) are all mutually consistent with one another. It is 
important to note that we do not conclude that gradients in the 
vicinity of the earth are in fact large at low energies, but only 
that on the basis of the model used large gradients and small aniso- 
tropies are not inconsistent with each other. 

Since there is some evidence that low energy protons (T -w 20 MeV) 
might be of solar origin (e.g. Kinsey, 1967), which is not consistent 
with the above result, we point out thac there are other possibilities. 
For example, (i) the model we have used might be inadequate, (ii)ft 
might not be large in the energy range in question, and (iii) the 
observations we have used might, for some reason, be incorrect. It 
is evident that to establish that the low energy particles are indeed 
of solar origin is not an easy matter, and if we are to use arguments 
similar to those developed in this paper, observations of the local 
spectrum, anisotropies, and radial gradients should be performed simul- 
taneously and with considerable care. 
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4 , Discussion 

Previous studies of the relationship between the radial gradient 
and radial anisotropy of cosmic rays have been carried out by Gleeson 
and Axford (1968b), Jokipii and Parker (1968), and Forman (1968b), 

In each case an attempt was made to deal directly with Equations (1) 
and (2), or with equivalent forms, rather than with asymptotically 
valid aoproximations as we have done in the present paper. 

The first such study to be carried out was that of Gleeson and 
Axford (1968b), who used Equation (2) together with observations of 
the local spectrum and radial gradient of cosmic rays, and also with 
an estimate for K based on the spectrum of interplanetary magnetic 
field fluctuations (e.g. Jokipii, 1966), to predict the behavior of 
the radial anisotropy as a function of particle kinetic energy. They 
found that the radial anisotropy should be small and positive for 
kinetic energies greater than about 600 MeV/nucleon, negative in the 
range 40 - 600 MeV/nucleon, and positive for low energies. In fact 
this prediction is qualitatively in agreement with observations 
(McCracken and Rao, 1966; Rao, et al., 1967). Since the argument 
is based on an exact, not an approximate, equation involving quan- 
tities which can be determined locally, it should be possible in 
principle to make correct predictions provided the equation itself 
is a valid representation of the behavior of the cosmic rays. 

For high energies Gleeson and Axford (1968b) used essentially 
the same arguments as we have used in Section 2. For intermediate 
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and low energies they had to estimate the term K dU/dr using the 
gradients determined by O'Gallagher (1967), assuming K a 8P through- 
out the entire energy range and matching K to the high energy value 
obtained £rom Equation (5). The resulting value of K dU/dr is small 
compared with CVU (c.f. case (b) of Section 3) and hence the general 
nature of the predicted anisotropy is rather insensitive to errors in 
K, and dU/dr unless K is actually very much larger than estimated. The 
predicted magnitude of the anisotropy does not however agree with that 
found by Rao, et al. (1967) for the energy range (7.5 - 45 MeV.) . 

Jokipii and Parker (1968) used Equation (1) in an Integral form 
in their discussion of this problem. If Equation (1) is multiplied 
by f and integrated with respect to r successively by parts in the 
range (0, r), one obtains: 
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on putting S * 3vU| , and using the condition r*S ” 0 as r 0 which 
corresponds to having no sources or sinks at the sun (i.e. no solar 
cosmic rays). Jokipii and Parker make the assumption that the gradient 
in the vicinity of the earth can be extrapolated all the way back to 
the sun, so that the higher derivatives vanish and the anisotropy is 
given by: 
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It Is clear, however, that such an assumption cannot be valid if the 
gradient in the vicinity of r ■ 1 A.U. exceeds 1007. A.U. as reported 
by O'Gailagher (1967). It is obvious from simple order of magnitude 
arguments and from consideration of analytic solutions (Fisk and 
Axford , 1969) that in situations where the radial gradient is not 
small, the terms in the series on the right side of Equation (18) 
becoming increasingly large in magnitude as i increases. Accordingly, 
one cannot expect Equation (19) to correctly predict the radial aniso- 
tropy using the large radial gradients reported by O'Gallagher (J967), 
and one cannot assert on the basis of such a prediction that there is 
any inconsistency between these radial gradients and the radial aniso- 
tropies observed by Rao, et al. (1967). 

At high energies, where the radial gradients can be expected to be 
small, Equation (19) is a useful result and indeed is the same as 
Equation (7) to within a small numerical factor. It is evident from 
consideration of analytic solutions (Fisk and Axford, 1969) that this 
numerical uncertainty arises because the terms r^d^U/dr* are in general 
of comparable magnitude when ft < 1 and hence one depends on the coef- 
ficients (•!)*•/ (1+2) ! for convergence of the series in Equation (18). 

In some simple cases the series can be summed (Fisk and Axford, un- 
published), and in particular if K *Ko r » Equations (7) and (18) 

(but not (19)) yield the same result for ft «1. 

Forman (1968b) has approached the problem from a somewhat different 
point of view, and points out that if the radial gradient of the intensity 
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is large at low energies as suggested by the observations of O'Gallagher 
(1967), then the radial gradient of the anisotropy, (r/5) (dg/Br) , must 
also be large at these energies. This result can be deduced easily 
from the approximate equations derived in Section 3. Thus, unless 
there was a chance cancellation, one might have expected a noticeable 
variation in $ even within the small variation of heliocentric distance 
(0.8 * 1.1 A.U.) involved in the observations of Rao, et al. (1967), 

On this basis the predictions of the model and the observations of 
the anisotropy and of the radial gradient are mutually inconsistent. 
There is no contradiction with the results given in Section 3 of this 
paper, however, since the comparisons involve different quantities. 

In effect, Forman's argument confirms our use of Equation (9); our 
neglect of 2S/r compared with dS/dr is essentially equivalent to 
assuming that dg /dr » §/ 
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FIGURE CAPTIONS 

Fig. I A plot of the radial gradient of the intensity determined 
by Equation (5) and the radial anisotropy determined by 
Equation (7), using the spectrum for protons given by 
Gloeckler and Jckipii (1967). K is taken to be propor- 
tional to 3P and is chosen to have a value K ■ 9.3 x 1C? 1 
cm 3 sec’ 1 for protons with T ~ 5 GeV in order that the 
predicted gradient at this energy is ~ 9%, in agreement 
with the observations of O'Gallagher (1967). V is taken 
to be 400 km sec* A . The Compton-Getting factor C cor- 
responding to this spectrum is shown by Gleeson and Axford 
(1968b). 

Fig. 2 A spectrum which provides a reasonable fit to the available 
observations of low energy proton intensities in 1966. The 
symbol 0 is used to represent the observations of Fan, et al. 
(1968), and the symbol A, the observations of Badhwar, et al. 
(1968). In case (c) discussed in the text, the spectrum is 
assured to result primarily from solar protons at energies 
below 10 MeV, and primarily from galactic protons at ener- 
gies above 60 MeV; the two dashed curves shown in this figure 
represent solar and galactic proton spectra, into which the 
total spectrum can be decomposed in the (10 - 60 MeV) energy 
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Fig, 3 A plot of the radial anisotropy determined by Equation (10) 
or (13) for each of the three cases discussed in the text, 
vs T. 

Fig, 4 A plot of the quantity T) determined by Equation (12) for 
case (a) discussed in the text, vs T. 

Fig. 3 A plot of the radial gradient of the intensity determined 
by Equation (16) for case (b) discussed in the text, vs T. 
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ABSTRACT 

Approximate equations which describe the behavior of cosmic rays 
in the interplanetary medium under suitable conditions are used to make 
comparisons between observations and theoretical predictions of radial 
gradients and radial anisotropies. In the high energy region there appear 
to be no inconsistencies between theory and observations. In the low 
energy region it is shown that theoretical predictions of the radial 
anisotropy expected from jlarge radial gradients of the intensity are not 
inconsistent with observed radial anisotropies. However, in the latter 

J ' 

case there are other inconsistencies, which suggest that some Aspects 

i 

of the observations or of the theory (or both) are unsatisfactory. 

1 . Introduction 

The modulation of cosmic ray particles in the interplanetary 
medium can be discussed in terms of the quasi-steady, spherically- 
symmetric model developed by Parker (1965, 1966), Gleeson and Axford 
(1967, 1968a), and Jokioii and Parker (1967). In this paper observa- 
tions of radial gradients and radial anisotropies of both high energy 
and low energy cosmic rays are compared with the predictions of the 
approximate equations derived for this model by Gleeson and Axford 
(1968b, c) and by Fisk and Axford (1969). 

It has been shown by Gleeson and Axford (1967) that the cosmic ray 
number density U(r,T) and streaming (or radial current density) S(r,T), 

> ■ 

per unit interval of kinetic energy T, satisfy the equations: 


